clear all 
clc
m=[0.2973    0.3964    0.1712    0.1351
    0.2703    0.4144    0.1622    0.1532
    0.2703    0.4505    0.2162    0.0631
    0.4234    0.1802    0.1081    0.2883
    0.2342    0.4234    0.2342    0.1081
    0.3514    0.3964    0.1261    0.1261
    0.3514    0.3604    0.0991    0.1892
    0.2793    0.3694    0.1622    0.1892
    0.2072    0.4324    0.2072    0.1532
    0.1818    0.4091    0.2727    0.1364];
a=[ 0.3545    0.1000    0.0455    0.5000
     0.3273    0.1455    0.0273    0.5000
     0.2545    0.1273    0.1000    0.5182
     0.3000    0.1182    0.0818    0.5000
     0.2909    0.0636         0    0.6455
     0.3636    0.0909    0.0818    0.4636
     0.3545    0.1364    0.2455    0.2636
     0.2909    0.0909    0.1182    0.5000
     0.2182    0.0727    0.1455    0.5636
     0.2000    0.0636    0.1727    0.5636 ];
 for i=1:10
   b(i,1)=fun3(a(i,1),a(i,2),a(i,3),a(i,4));
   c(i,1)=fun3(m(i,1),m(i,2),m(i,3),m(i,4));
 end
c 
b
